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Abstract 

We discuss the potentials of several experimental configurations dedicated to direct 
measurements of charged cosmic ray (CR) nuclei at energies > 100 GeV/n. Within 
a two-zone propagation model for stable CRs, we calculate light primary and sec- 
ondary nuclei fluxes for different diffusion and acceleration schemes. We show that 
the new detectors exploiting the long and ultra long duration balloon flights could 
determine the diffusion coefficient power index 5 through the measurement of the 
boron-to-carbon ratio with an uncertainty of about 10-15 %, if systematic errors are 
low enough. Only space-based or satellite detectors will be able to determine 5 with 
very high accuracy even in the case of important systematic errors, thanks to the 
higher energy reach and the less severe limitations in the exposure. We show that 
no uncertainties other than those on 5 affect the determination of the acceleration 
slope a, so that measures of light primary nuclei, such as the carbon one, performed 
with the same experiments, will provide valuable information on the acceleration 
mechanisms. 
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1 Introduction 



Charged particles arriving at Earth with energies between about 10 7 eV and 10 15 
eV are believed to have galactic origin. Their acceleration is very likely due to the 
action of supernova remnants, while their subsequent diffusion in the Galaxy is 
driven by the turbulent, irregular component of the galactic magnetic field. The 
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most abundant galactic cosmic ray particles have energies in the 100 MeV/n - 10 
GeV/n range, where several additional phenomena such as electromagnetic energy 
losses, convection, diffusive reacceleration and solar modulation are believed to con- 
tribute in shaping their spectra. The most realistic description of CR propagation is 
given by diffusion models, in which the several free parameters inherent to a specific 
model need to be fixed by observations. 

The low-energy tail of the CRs spectrum is shaped by several competing effects, so 
that it is very difficult to disentangle each physical component. On the other hand, 
most data have been collected for energies < 50 GeV/n and their interpretation has 
not yet led to a clear understanding of any of the above-mentioned physical ingredi- 
ents. The high energy part of the spectrum - let's say > 10 2 GeV/n - is basically due 
to acceleration and diffusion, all the other effects being minor if not negligible. The 
best-founded assumptions on the diffusion coefficient and the injection spectrum 
are power-laws, and this is what is grossly observed for the flux of nuclei. 
A wealth of experimental measurements are available, with different degrees of ac- 
curacy in the region up to 100 GeV/n; on the contrary, the higher energy region is 
poorly known. The most recent direct measurements in this region have been pro- 
vided by the series of balloon flights of JACEE [1] and RUNJOB [2], while most of 
the data results from indirect measurements, by means of ground arrays observing 
air showers. 

New data have recently been added by the ATIC Collaboration [3], connecting the 
lowest energy region to the highest energy available data. The new CREAM project 
[4] has been developed and is now in data taking phase, with the aim of dramatically 
increasing the available statistics in the energy region up to 500 TeV and possibly 
more, thus making new and precise measurements of the spectral characteristics of 
CRs. 

In the present paper, we explore the performances required for new detectors to 
disentangle the fundamental parameters describing the propagation of nuclei at en- 
ergies above about 100 GeV/n and below the knee region. In Sect. 2 we report 
on the status of experimental projects aimed at the measurement of CR nuclei at 
energies > 1 GeV/n. In Sect. 3 we describe a diffusion model for galactic CRs, and 
highlight the main features affecting the propagation of particles in the higher en- 
ergetic range. Sect. 4 summarizes the method employed for the simulation of the 
experimental conditions. In Sect. 5 we present our results focusing on the possi- 
bilities to measure the diffusion coefficient slope by means of the boron-to-carbon 
(B/C) data and the acceleration spectrum by means of primary light nuclei fluxes 
such as the carbon one. In Sect. 6 we draw our conclusions. 



2 Experimental status and projects 

In Table 1, we collected the past and current experiments as well as new projects 
directly measuring charged CRs nuclei in the energy range above 1 GeV/n and 
spanning a wide range in Z. 

Direct measurements of charge and energy of CR nuclei started in the sixties, ex- 
ploiting the balloon and spacecraft techniques. Scintillators and Cerenkov detectors 
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Table 1 

Balloon and space instruments and new projects. Experiments covering an energy 
range below 1 GeV/n are not considered in the Table. For all the new projects on 
satellite we assume 1000 days of data recording. At least 10 events are expected to 
be detected by ATIC, CREAM and ACCESS above the maximum quoted energy. 



Experiment 


Years 


Elements 


Energy 


r 


Exposure 








GeV 


m 2 sr 


(m 2 sr days) 


Balloon instruments 












MUBEE [12] 


1975-87 


1 < Z < 26 


10 4 -3 • 10 5 


0.6 


22 


SANRIKU [131 


1987-91 


8 < Z < 26 


2 -10 2 -3- 10 4 


4.81 


4.45 


JACEE [11 


1979-95 


1 < Z < 26 


2 -10 3 -8- 10 5 


2-5 


65-107 


RUNJOB [2] 


1995-99 


1 < Z < 26 


10 4 -3 -10 5 


1.6 


43 


TRACER [14] 


2002-04 


5 < Z < 28 


1.6 -10 2 -1.6 -10 4 1 


5 


70 


ATIC [3] 


2000-03 


1 < Z < 28 


10-10 5 


0.25 


3.5 


CREAM [4] 


2004^ 


1 < Z < 28 


10 3 -10 6 


0.5/1.3 2 


35-140 3 


Space instruments 












PROTON 1-4 [101 

L J 


1965-68 


All particles,p,He 


100-10 6 


0.05-10 


5-2000 


SOKOL [16] 


1984-86 


1 < Z < 26 


2 -10 3 -10 5 


0.026 


0.4 


HEAO-3 [11] 


1979-80 


4 < Z < 28 


9.6-560 1 


0.14 


33 


CRN [17,18] 


1985 


4 < Z < 26 


7- 10 2 -3 TO 4 


0.1-0.5/0.5-0.9 2 


0.3-3 


Projects 












NUCLEON-KLEM [19] 




1 < Z < 26 


10 2 -10 7 


0.19 


190 


AMS-02 [26] 




1 < Z < 26 


1-10 4 


0.5 


500 


ACCESS [23] 




1 < Z < 28 


10 3 -5 -10 6 


1/8 2 


300-8000 


PROTON 5 [20] 




1 < Z < 28 


10 3 -10 7 


18 


18000 


INCA-CSTRD [21] 




1 < Z < 28 


10 5 -10 7 


48 


48000 



(1) the energy range refers to oxygen; 

(2) the two numbers refer to low and high Z respectively; 

(3) CREAM had a succesful 40 days flight in December 2004-January 2005, but it 
is planned to exploit ULDB flights of ~ 100 days in the near future. 



were employed for charge measurement, while ionization calorimeters measured the 
nucleus energy [5,6,7,8,9,10,11]. 

Four experiments [1,2,12,13] exploited passive detectors like nuclear emulsions and 
X-ray films: due to their sensitivity also to inclined tracks, they allowed to reach 
a much bigger exposure compared to active detectors, thus in principle exploring 
the highest energy region (about 10 to 1000 TeV). The most precise measurements 
on relative abundances of elements with Z from 4 up to 28 came from the HEAO-3 
mission, whose data have been extensively used to tune and check models. AMS-01 
(a spectrometer flown on the space shuttle Discovery including a permanent mag- 
net, a tracker, time-of flight hodoscopes and a Cerenkov counter) brought further 
information on proton and helium up to about 100 GeV/n [15]. 
Among the present experiments, a fully active bismuth germanate (BGO) calorime- 
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ter flew on ATIC [3], while TRACER [14] employed a transition radiation detector. 
The development of the Long Duration technology is now giving the possibility to 
balloon experiments to fly for more than 20-30 days, thus reaching much higher 
exposures. CREAM [4] employs for the first time both a calorimeter and a transi- 
tion radiation detector, thus allowing an intercalibration of the energy measurement 
with two independent techniques with different systematic biases. This apparatus is 
planned to exploit the Ultra Long Duration Balloon flights technology, which relies 
on new superpressure closed balloons with a pumpkin-shape design and is now in 
test phase; it will hopefully fly up to 100 days. 

New space instruments have also been proposed: a) the NUCLEON russian satellite 
program which will include the KLEM detector (silicon microstrips and scintillator 
strips) [19]; b) a ionization calorimeter on a series of Proton-5 satellite missions [20], 
with the aim of measuring the spectrum of elements up to more than 1000 TeV, 
reaching an accuracy of 2-5% in the slope determination; c) INC A [21], a ionization 
neutron calorimeter with 48 m 2 sr acceptance; its main aim being the study of the 
electron spectrum, it is expected to measure also primary nuclei in the knee region. 
A proposal for an integration of this calorimeter with a Compton scatter transition 
radiation detector has been recently presented [22]. The ACCESS [23] project was 
not chosen for the last MIDEX phase by NASA, but has been included in the hope 
of a reproposition as a free-flyer. 

We should also mention the sophisticated spectrometers which were flown on bal- 
loons, with the main aim of measuring antimatter and light isotopes [24]. They are 
not included in the Table because of their quite different aim, but they also gave 
very precise information on the proton and helium flux up to around 100 GeV/n; 
the energy range has been extended by the Bess- TeV [25] upgrade of Bess-98, reach- 
ing about 500 GeV and 200 GeV/n for proton and helium nuclei respectively. The 
new AMS-02 [26] large acceptance magnetic spectrometer has been conceived to 
study origin and structure of the dark matter and to measure antinuclei; it will also 
be able to extend the knowledge on the composition of charged CRs from the 100 
MeV/n region up to about 1 TeV/n. 



3 Propagation of CRs in the Galaxy 

Despite scarce theoretical and observational knowledge of the key parameters re- 
sponsible for acceleration and propagation of galactic CRs, phenomenological mod- 
els able to reproduce data can be built. The most realistic propagation models are 
the diffusion ones, even if the so-called leaky box model has been often preferred 
in the past for its simplicity [27]. In Refs. [28,29], a two-zones diffusion model has 
been developed and shown to reproduce several observed species in the low-energy 
part of the CR spectrum (~ 0.1-100 GeV/n). In this model, the Galaxy has been 
cylindrically shaped (with R=20 kpc), with a thin disc (half-hight h=0.l kpc) con- 
taining the sources and the interstellar medium (ISM) surrounded by a diffusive 
halo of half-thickness L ~ 2-15 kpc. The transport equation for the nucleus j in 
a diffusion model is in principle valid in a very wide range of energies and can be 
written as: 
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(1) 



+2h6(z) qlQ ] {E)q{r)+Y J ^ k3 N k {E,r^)-T^N\E,r^) 



= 2h5(z) — V(E)Ni(E,r,0) - d\E)— N*(E, r, 0) 




fc=i 




Steady-state has been assumed and N 3 (E, r, z) is the differential density of the nu- 
cleus j as a function of energy E and galactic coordinates (r, z). In this equation, 
the first term represents spatial diffusion with diffusion coefficient K(E), which 
has been assumed to be independent of Galactic coordinates. V c is the convection 
velocity, assumed here to be constant throughout the Galaxy (except in the thin 
disk) and to be directed outward along the z-direction. The term proportional to 
2h in the left-hand side of Eq. (1) takes into account all the sources of cosmic rays: 
primary sources with injection spectrum Qi(E), secondary spallative production 
from heavier nuclei and destructive reactions (radioactive-decay terms have been 
omitted for simplicity). The right-hand side of Eq. (1) contains - through the coeffi- 
cients W(E) and d^(E) - the terms responsible of the energetic changes suffered by 
charged particles during propagation: Coulomb, ionization and adiabatic expansion 
losses, and gains due to second order reacceleration. Diffusive, or continuous, reac- 
celeration is due to the scattering of charged particles on the magnetic turbulence 
in the interstellar hydrodynamical plasma. The diffusive reacceleration coefficient 
is related to the velocity of such disturbances, called the Alfven velocity Va, and is 
naturally connected to the space diffusion coefficient K{E). For example, diffusive 
reacceleration contributes significantly in shaping the boron-to-carbon (B/C) ra- 
tio at kinetic energy per nucleon E around 1 GeV/n. Indeed, all energy losses and 
gains are effective only in the low-energy tail of the CR spectrum and thus are irrel- 
evant in the analysis carried in the following of our paper, which deals with E> 100 
GeV/n. The so-called sporadic, or distributed, reacceleration has been considered 
in the literature [30,31] and is based on the possibility that CRs be reaccelerated 
during their wandering in the Galaxy by supernova remnants (SNRs) and gain a 
small amount of energy at each encounter. 

Convection may dominate at low energies, depending on the value of V c and if S is 
large. It may compete with diffusion up to few tens of GeV/n, but its role becomes 
negligible at higher energies. The importance of spallation processes depends on 
the nucleus, but at energies > 100 GeV/n they are not much relevant except for 
heavy nuclei (e.g. iron) [32]. We also note that the different assumptions on the 
distribution of the ISM - homogeneously distributed (as assumed in our analysis) 
or spatial dependent - are irrelevant at high energies and for the species considered 
in our analysis. 

Diffusion depends on the rigidity R (R = p/Z) of the particle and the diffusion 
coefficient is usually assumed to have the form: K(E) = K~o[3R & . K$ is linked to 
the level of the hydromagnetic turbulence and 5 to the density spectrum of these 
irregularities at different wavelength. The Kolmogorov theory for the turbulence 
spectrum predicts 5=1/3, while the case for a hydromagnetic spectrum with 5=1/2 
has been obtained by Kraichnan [33]. On the other side, the phenomenological inter- 
pretations of CR spectra have not led to fix 5 because of the complicated treatment 
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of the physical phenomena which are relevant at the energies in which most cosmic 
particles have been collected. Fits to the B/C ratio within the diffusion model of 
Eq. (1) prefer high values for 6 (~ 0.5-0.75), while the Kolmogorov spectrum turns 
out to be disfavored [28,29,34]. The 5=1/3 case can reproduce quite well the pe- 
culiar B/C peak observed at ~ 1 GeV/n but it tends to overestimate the data at 
increasing energies, predicting a quite flat B/C ratio (see also later). The highest 
values of 5 seem on the other side not very realistic. In fact they are only marginally 
consistent with observations of interstellar scintillation [35,36] and are highly incom- 
patible with the level of anisotropy measured at 1-100 TeV [37]. Data on CR fluxes 
do not extend to sufficiently high energies and are not sufficiently precise in order 
to constrain 5. 

The acceleration spectrum of primary CRs is believed to be determined by super- 
novae (SN) remnants and super-bubbles [38], which are the only known engines in 
the Galaxy able to provide the right amount of energy to particles up to ~ Z x 10 14 
eV. The acceleration spectrum follows a power-law in momentum, Q(E) oc p~ a , 
with a located somehow between 2.0 and 2.5. Even if a precise value for a cannot 
be predicted, the effective spectral index as derived by several indirect observational 
tests and by numerical studies is close to a ~ 2.0-2.1 [38] (and refs. therein), [39]. 
Very recently, the TeV 7-ray image of a SN remnant - already observed in the X-ray 
spectrum with a very similar morphology - strongly indicates that high energy nu- 
clei are accelerated in this site [40] . The 7-ray data are well reproduced by a photon 
spectral index which seems to prefer low a values, even if the experimental sensi- 
tivity is not yet sufficient to put strong constraints. Higher values of a ~ 2.4 are 
favored by propagation models for GCRs with very low S, such as the Kolmogorov 
spectrum. 

The above-described diffusion model has been tested on the B/C and sub-Fe/Fe 
ratios and shown to fit well existing data, which we remind lie in the low-energy 
tail of the galactic CR spectrum [28,29]. The propagation parameters able to repro- 
duce such ratios give also rise to a secondary antiproton flux which is in very good 
agreement with measurements [41]. A further, yet weaker, degree of consistency is 
given by an analysis of radioactive isotopes [42]. In this case, the diffusion model 
has been modified in order to take into account a local under-dense region, rele- 
vant when dealing with short-living species. The model has also been validated at 
higher energies through a study of the mean logarithmic mass of the CR [32] beam. 
For further details on this propagation model and considerations on other possible 
approaches we refer to Ref. [43] and references therein. 



4 Detector design and simulation 

To pursue the objective of extracting information about crucial propagation param- 
eters in the high energy region, an experimental apparatus should have the following 
characteristics: 

(i) A large exposure, to have enough statistics: to detect about 10 nuclei at energies 
above 10 TeV/n, if taking the elemental fluxes from [44], a collecting power T ~53 
m 2 sr days is to be reached for carbon nuclei, while more than 250 m 2 sr days are 
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necessary for Iron. 

(ii) A good energy resolution: this is especially crucial to detect changes of slope 
in the energy spectra (e.g. the knee in the proton spectrum) or to study spectral 
smoothness. A constant 40% energy resolution is enough to see an increase in the 
spectral slope of about 0.3 crossing the knee. The same resolution or better is re- 
quired to find small deviations from smoothness, of the order of 10% [20]. 
(hi) A charge resolution such as to distinguish e.g. carbon from boron nuclei; models 
ask for a B/C ratio of some percent around 1 TeV/n, so that a resolution of about 
0.2 charge units is needed. 

Combining different detectors in the same experiment gives a powerful tool to over- 
come individual technical limitations: for example, redundant measurements of en- 
ergy allow a cross-calibration of the detectors, thus overcoming the problem of direct 
energy calibration which at the highest energies explored is possible for TRDs but 
not for calorimeters. 

A very careful study of systematic errors is mandatory, as these are the main cause 
of uncertainty in all the measurements. Exposures and efficiencies in selecting and 
tracking the events can be studied by means of Monte Carlo simulations; for both 
energy and charge determination, redundancy can help in reducing the ambiguities. 
Although a serious discussion of systematics depends on each specific detector, a 
constant contribution to the overall error on the flux measurements will be intro- 
duced in our analysis. 

In order to explore the performances of new detectors, a simulation has been built 
with the following ingredients: 

(a) the input flux; the particle spectrum is numerically given by the propagation 
model. Different power laws with increasing slope in contiguous energy ranges are 
fitted to the spectrum to better follow the expected behavior. The flux is normalized 
to the best experimental low energy datum (e.g. the flux as derived by [11] at 10.6 



(b) the experimental apparatus, in terms of collecting power and energy resolution. 
We use r = 1.3, 5 and 10 m 2 sr, similar to the quoted geometrical factors of some of 
the current and future experiments and a set of exposure times of 30, 100 and 1000 
days, which roughly correspond to long and ultra long duration balloon flights and 
satellite conditions, respectively. 

The expected number of events above a given energy E for an assumed collect- 
ing power T and input spectrum dN/dE = CqE~^ m.~ 2 s^ 1 sv^ 1 (GeV /n)^ 1 can be 
written as: 



where C\ depends on the slope and normalization of the input flux. 

Following the approach of [45], the cumulative distribution function for E in the 

considered energy range can be written as 



GeV/n). 



N(> E) = rCi£~ 7+1 



(2) 



&(E) = 1 - 



N(> E) - N{> E b ) 
N(> E a ) - N(> E b ) 



[E a <E< E b ] 



(3) 
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and the probability distribution function for the events is obtained by differentiation 



The number of events N of Eq.2 is the mean value of the Poisson distribution of 
the true number of events in a given range [E a — Ef,]; the poissonian fluctuation 
is computed in small energy bins (equal bins in logarithmic scale so that, for any 
interval j, Ej + \/ ^=constant) covering the full interval, in order to correctly weight 
also the bins with few events. The energy of each i-th event is then randomly sampled 
in the energy interval [E a — Ef,] from the power law spectrum; being a standard 
uniform variate 

Ei 

$(£.) = J d<S>(E')/dE'dE' = n (5) 

E a 



and 



Ei = {nm 



-7+1 



e: 



-7+1 



)+E~ 



7+1}- 



(6) 



In each interval [Ej-Ej+i], the mean energy is 



7 u h j ~ h 3+i 



<E>=(^)( ' '7 ) ( 7 ) 



The detector response is modeled as a Gaussian distribution with mean equal to the 
input energy and width given by the energy resolution. In the following, a constant 
40% energy resolution will be assumed; the case for resolutions either decreasing 
or increasing with energy will however be checked. The detection efficiency is here 
assumed not to depend on energy and will be included in the systematic uncertain- 
ties. 



5 Results and discussion 



The aim of this section is to study the potentials of different experimental configu- 
rations in connection with the determination of the diffusion coefficient slope 5 and 
the acceleration power spectrum a. We show how well their physical value could be 
singled out with high-energy data on the B/C ratio and the carbon flux. 
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Fig. 1. B/C ratio as a function of kinetic energy per nucleon. Solid curves: theoretical 
predictions corresponding, from top to bottom, to cases 1-2-3-4-5, respectively (see 
text for details). Dashed line: simple power law spectrum B/C=1.4 E~°' eo . 

5.1 B/C and the diffusion coefficient slope 



The boron-to-carbon ratio has always been considered the best quantity to study 
diffusion properties. At energies around GeV/n, this ratio is strongly affected by 
the low-energy phenomena described in Sect. 3. In particular, the observed bump at 
~ 1 GeV/n can be naturally explained by diffusive reacceleration, without invoking 
any artificial break in the diffusion coefficient. 

The most precise B/C data have been obtained in the energy range 0.6-35 GeV/n 
by the HEAO-3 experiment [11], the highest energy points having been collected by 
CRN [17]. In Fig. 1 we plot a collection of data on the B/C ratio giving emphasis 
to the higher energy part of the spectrum. Along with the experimental data, we 
also plot the theoretical predictions calculated with the full diffusion model briefly 
outlined in Sect. 3. 

We have chosen five illustrative configurations of the propagation parameter space, 
which will be also employed in the following of our analysis as bench marks. At 
the energies we deal with, the parameter which mostly shapes the B/C ratio is the 
diffusion coefficient slope, namely 5. The constant K$ in the diffusion coefficient 
enters with the source abundances in the global normalisation of B/C, so that its 
precise value is of scarce relevance in the B/C predictions. We thus identify cases 
1-2-3-4-5 with 5 = 0.3, 0.46, 0.6, 0.7, 0.85, respectively. All the other propagation 
parameters (Kq, L, V c , Va, a) have been fixed according to the best fits to B/C data 
[29] . Indeed their exact value is not very important when working with a secondary- 
to-primary flux ratio at energies > 100 GeV/n. The dashed line in Fig.l corresponds 
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to the power law 1.4E7 , where the coefficient 1.4 has been fixed to the theoret- 
ical prediction for the 5=0.6 case at 1000 GeV/n. This line shows that for energies 
> 1000 GeV/n the B/C ratio can be well approximated by a power law in en- 
ergy with index —S, as one expects from the high energy limit: secondary /primary 
oc K(E)^ 1 oc R~ s . At E=100 GeV/n the difference between the simple power law 
and the calculated flux is < 20%. Below this energy, it is clear that the effects of 
energy changing, convection, halo size and spallations become more and more im- 
portant in shaping the B/C ratio, as studied in Refs. [28,29]. 

This figure wants to emphasize the discrepancy among the models - basically due 
to the 5 values - with increasing energy. At energies ~ 100 TeV/n the predictions 
for B/C ratio in the 5=0.3 and (5=0.85 case - both not excluded by low energy data, 
even if (5=0.85 is a quite extreme value - differ by two orders of magnitude. Different 
predictions for 5 lead to B/C whose relative values increase with increasing energy. 
This high energy region is thus very interesting from an experimental point of view, 
since it could bring - at least in principle - to a clear determination of the actual 
diffusive regime. 

Finally, we want to underline that the injection spectrum has no relevance in the 
calculation of this ratio. We have checked that for E > 10 GeV/n the results for 
the B/C ratio are practically unchanged if we vary a (equal for each nucleus) in 
the very large range 1.8-2.5, as well as if we fix a for each nucleus (a- 7 , indeed) as 
derived in Ref. [44]. 

Following the procedure outlined in Sect. 4, the expected number of carbon and 
boron events can be derived by selecting a given set of experimental conditions 
(exposure, time and resolution) and an input flux. Fig. 2 shows the B/C ratio as ob- 
tained from case 1 (upper panel, corresponding to (5=0.3) and from case 3 (bottom 
panel, (5=0.6), for three different collecting powers. The maximum energy at which 
the B/C ratio can be measured with a significant number of events - we are requir- 
ing at least 10 events for Boron nuclei - goes from about 900 GeV/n to more than 
10 TeV/n when the exposure increases from about 40 m 2 sr days up to 10000 m 2 sr 
days. The various correction factors which must be considered in the analysis of real 
data (e.g. due to selection efficiency or interaction losses in the apparatus) would 
have the effect of lowering the maximum detectable energy. As an example, a global 
efficiency of 30% would shift it by a factor 1.5 for the smallest exposure considered 
here and case 3 as input. The same effect is produced for balloon experiments due 
to spallation processes in the residual atmospheric grammage of 4-5 g/cm 2 . 

/,From this figure it is evident that the two theoretical models could be easily dis- 
tinguished already by experiments with acceptance of about 1 m 2 sr and flying time 
of forty days (in our example: 1.3 m 2 sr and 30 days data taking). Note however 
that the figure has been obtained without considering the systematic uncertainties. 
The inclusion of a finite energy resolution of the detector in the simulation has 
important consequences on the ability to measure the B/C ratio and to distinguish 
among different models: it produces in fact a distortion which reflects the change 
in the carbon and (more) in the boron spectra due to the energy fluctuation. The 
steepness of the spectra makes it more likely to move low energy events to high 
energy bins, so that the effect is bigger for weaker spectra. In the case of model 1, 
for example, the measured B/C will be some percent higher than what expected for 
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Fig. 2. B/C ratio as expected from input case 1 (upper panel) and 3 (lower panel) 
for different experimental exposures. The expected ratios from cases 1 to 5 (from top 
to bottom, respectively) are also shown. 



an ideal detector, that is a "zero resolution" one, almost constantly from 100 GeV/n 
on; on the contrary, using the input spectrum of model 5 the B/C ratio appears to 
be from 12 to 18% higher when going from 100 to 1000 GeV/n. A correction for 
this effect has to be applied before the comparison with expectations; the binning 
of the data should be carefully chosen in order to be bigger than the uncertainty 
due to energy resolution. 

Using two different input models j and k, we simulated the expected number of 
carbon and boron events. The difference between the B/C ratios so derived can be 
quantified by performing a relative x 2 test 



X 



2 (( B / C )ij ~ (B/C) ik j 

^r— ^ (8) 
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where we put at denominator the variance of the difference. In order to fully exploit 
the high energy region which can be explored by the current and future apparata, 
where only acceleration and diffusion dominate the cosmic rays spectrum, the sum- 
mation starts from 100 GeV/n. 

The ratio between the carbon and boron fluxes is only mildly sensitive to system- 
atic errors (most of which are supposed to be independent from the considered 
nucleus), which are however to be included in the total error. Since a full discussion 
on systematics can only be done for each specific apparatus, we consider here only 
a constant systematic error to be added in quadrature to the statistical one. 
Fig. 3 shows the result for the data sets respectively expected from model 3 and 
each of the other ones as a function of the systematic uncertainty, for three different 
experimental collecting powers. We should remind that even at the lowest energy of 
100 GeV/n the percentage difference between models is big, going from about 7% for 
models 3 and 4 up to 30% for models 1 and 2. Exposures are such that even for the 
smallest experimental configuration considered here the statistical errors are quite 
small and in the ideal case of no systematic effects all models could be separated. 
As obvious, the sensitivity increases with differences between the two compared 5 
values, in our case the ones labeled (1-3) and (3-5) for which (Sj — <5fc)=0.30 and 
0.25 respectively. For example, a 5=0.3 case can be distinguished from a 5=0.6 case 
with 90% C.L. if the systematic error stays below 12% even in the less favourable 
experimental conditions considered. In the case of a satellite-like exposure, these two 
models could be distinguished even in presence of a 15% systematic error. On the 
other hand, the range of sensitivity raises up to more than 20 TeV/n for satellite- 
based experiments, thus allowing to tag bench marks with A<5=0.10-0.15 (e.g. case 
2 from 3 or case 3 from 4), but only if systematic uncertainties are kept below 10%. 

The analysis shown in Fig.3 is meant to illustrate the ability of a given experimental 
setup to discriminate between two fixed and arbitrarily chosen 5 values. In order to 
evaluate the accuracy in the determination of 5, we now perform a % 2 minimization 
procedure on the simulated B/C data, leaving 5, Kq, L and V c as free parameters. 
As expected, the resulting \ 2 has no structure with respect to Kq, L and V c , whose 
effect is in practice reabsorbed in a global normalization of the spectra. On the 
contrary, the calculated x 2 s are such that a minimum for 5 can always be clearly 
identified. 

If systematics can be considered negligible, the determination of the diffusion coef- 
ficient slope can be obtained at 10-15% significance from B/C data collected with 
the 39 m 2 sr days simulated experiment, and at 8%-10% for the 500 m 2 sr days 
case. However, we find that 5 can be determined with significance level of 15% or 
10% respectively for the 500 m 2 sr days or 10000 m 2 sr days simulated experiment, 
if the systematic errors are of order 10%. Very low accuracy is predicted in this case 
from the lowest exposure experiment (39 m 2 sr days). 

As expected, the obtained accuracies are consistent with the probability curves 
shown in Fig.3; e.g., the case 3-4 presented in this figure (bottom left) considers two 
5 differing by 15-20% and shows that they are perfectly separated only if the sys- 
tematic uncertainties lay below 2-3% for the smallest collecting power. This agrees 
with our previous conclusion, namely that the accuracy on 5 is at the level of 10% 
for a collecting power of 39 m 2 sr day and negligible systematics. 
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Fig. 3. Probability of distinguishing between two input models vs system- 
atic error for three different collecting powers. Remind that, for cases 1-5, 
8 =0. 30, 04 6, 0. 60, 0. 70, 0. 85 respectively. 

We emphasize once again that the above results are derived including systematic 
uncertainties which are assumed to be independent on energy: while this is true 
for trigger efficiency or event reconstruction, or even for Monte Carlo corrections, 
other systematic errors like those due to energy resolution will depend on energy 
and could alter the considered fluxes. 

A spectral deformation of the B/C ratio at high energies has been proposed in Ref. 
[31] as due to distributed reacceleration (see above). CR particles can have sporadic 
encounters with SNRs, during which they could gain small amount of energies. The 
modification of the primary spectrum comes out to be independent of energy, while 
the effect for secondaries increases with energy just as the diffusion coefficient in the 
Galaxy. Thus, the ratio of reaccelerated B/C re to standard B/C would be roughly 
described by (B/C) re / (B/C) cx (a + S - l)/(a + 26 - 1)(p/pgcr) S , where p GC R is 
the momentum corresponding to the maximum of the observed particle spectrum 
(in Ref. [31] it is assumed Egcr = 0-Q GeV/n). The real situation is rather more 
complicated than this description, which nevertheless gives us the correct trend and 
the possible magnitude of the effect. The expected deformation of the B/C ratio 
would be greater for greater 5. Moreover, as figured out in Ref. [31], its intensity 
would strongly depend on the value of the circumstellar hydrogen density, Njj- The 
higher Njj, the lower the distortion. In the case of 5=0.6, an effect of one order of 
magnitude would be expected at 1 TeV/n for iV#=0.003 cm" 3 , a very low number, 
indeed. For higher and more plausible values of the hydrogen density, the B/C ratio 
could be enhanced by a factor of two. If the diffusion is Kolmogorov-like, only very 
small Njj could give an observable effect (about a factor two). But in this case it 
would be less ambiguous to correctly interpret data, while in case of higher 5 this 
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sporadic reacceleration effect could be mimicked by smaller 5 effects. 
Non-stationary models have also been studied, taking into account the fact that 
SN might be considered discrete sources in the Galaxy, so that a time-dependent 
diffusion equation should be solved [49]. However, it seems plausible that the B/C 
ration would not be modified in shape, but only in a overall normalization factor 
[50]. 



5. 2 Primary fluxes and acceleration spectra 



Primary fluxes trace back to the diffusion and acceleration properties, since they 
are produced directly in the acceleration sites. Charged CR nuclei heavier than 
protons suffer catastrophic losses due to nuclear destruction on the interstellar H 
and He. These reactions become irrelevant on the CR spectrum at energies of about 
1 TeV/n, depending on the nuclear species. The effect of inelastic interactions lowers 
and flattens the flux, giving rise to a final shape that is harder than a pure power-law 
-E -7 , being 7 = a + 5. This effect has been studied by computing the carbon flux 
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Fig. 4. Ratio of calculated carbon fluxes from the same propagation model, where 
nuclear destructions have been turned off at the denominator. 

according to the full diffusion model (case 2 and a = 2.05, but in this context their 
value is not relevant) with and without inelastic destructions; the result is plotted in 
Fig.4 (for iron flux, see Fig. 1 in [32]). The difference between the two cases is about 
30% at 100 GeV/n, and at 1 TeV/n it is still 5%. Thus, it is in principle not correct 
to expect a pure power law spectrum from energies above few tens of GeV/n, as 
often claimed in literature. Of course, the error bars on the existing experimental 
data for primaries (except for protons, maybe) still justify this approximation. 
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^From the experimental point of view, data about fluxes of primary nuclei were 
collected in the past for Z<26 up to few TeV/n, but with decreasing statistical 
significance as the energy grows higher. At the highest energies, the only available 
data sets come from emulsion chamber experiments [1,2,16] which measure groups 
of elements; systematics are quite large, e.g. the results for the CNO group differ 
between JACEE and RUNJOB by a factor of 2 in normalization. More recently, 
the first data from a one day test flight of TRACER have been published for Z>8 
nuclei [48]. 

We have calculated the carbon flux for a propagation model in agreement with 
the B/C analysis presented above, with the only aim to qualitatively reproduce 
the data at the energies in discussion. So we do not care much about the level of 
agreement with lower energy data, aware that a thorough study of the propagation 
of primaries at low energies would require many deeper insights, which are beyond 
the scope of the present paper. Figure 5 shows a compilation of the data on carbon 
flux together with the calculated flux for model 3 (5=0.6) and a=1.9, 2.05 and 2.2. 
The normalization has an important role in this calculations and may be ascribed 
to the value of the source abundances, which are very poorly known, and of Kq. We 
refer our normalization to the fluxes measured at 35 GeV/n [11], where the total 
error (statistical and systematic) is about 15%. 
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Fig. 5. Data on carbon flux multiplied by E 2 ' 7 ; 3 theoretical models are shown, cor- 
responding to 5=0.6 and a=l. 9,2. 05,2. 2 from top to bottom, respectively. 

The expected number of events from the above described model and for different 
values of a have been calculated following the procedure described in Sect. 4 for 
different experimental exposures; the corresponding flux is shown in Fig. 6. 



We have also checked that this flux is almost independent of the assumed energy 
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Fig. 6. Expected C flux for input model with a = 2.05 and 5 = 0.6 and 3 different 
exposures. Theoretical curves as in Fig. 5. 



resolution by changing it from a constant value to an either decreasing or increasing 
behavior as a function of the particle energy (from 40% at 100 GeV/n to either 
50% or 30% at 1000 GeV/n). Indeed, it can be shown that the energy resolution 
plays a minor role as compared to the event statistics (to collecting power) in the 
estimation of spectral slopes for simple power laws [45]. 

The statistical error on the expected events is less than 2% up to about 500 GeV/n, 
rising to 10% at 2 TeV/n for the smallest exposure, with 500 m 2 sr days it lowers 
to few percent at 1-2 TeV/n: the new detectors will so have full ability to cover 
the high energy range with sufficient statistical precision. With a x 2 minimization 
procedure, the expected data are fitted to a power law and the slope 7 is determined 
with an error lower than 1% in a range AE=0.1-10 TeV/n. The effect of adding a 
constant systematic error of 10% is of more than doubling the error on 7. 
This result implies that the experimental errors on the carbon flux (the same con- 
clusion is reached for oxygen) will not add any uncertainty in the determination 
of the acceleration power law index a other than the one carried by the diffusion 
slope 5. If the latter is determined within 10% from B/C measurements, the same 
experiment will let us fix the former within an equal uncertainty level by means of 
primary fluxes determination. 

The knowledge of a will lead to a better understanding of the acceleration processes 
inside galactic supernova remnants. In particular, we note that even in a LDB mis- 
sion it could be possible to reach 1 TeV/n with quite good statistical accuracy, thus 
allowing to confirm or exclude whether a lies towards low- values 1.9-2.1, or in the 
higher range 2.3-2.4 at a 3 a level. 
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6 Conclusions 



A two-zone diffusion model for stable, galactic CRs has been extensively used to 
generate different predictions on the primary and secondary nuclei fluxes, which in 
turn are used to evaluate the expected events in various experimental configurations 
similar to present or next generation apparata. They will substantially exploit a high 
collecting power in the high energy region (> 100 GeV/n), where the low energy 
effects and most of the peculiarities of the model can be neglected. We have shown 
that new detectors exploiting the long and ultra long duration balloon flights, or 
even more future space missions, will provide valuable information about both the 
diffusion coefficient slope S and the acceleration spectrum a of galactic CRs. 
An experimental determination of 5 through the measurement of the B/C ratio is 
at reach of the new generation of detectors, some of which are already taking data. 
An experiment like CREAM - which during its first flight collected data for more 
than 40 days breaking all the previous flight and duration records 4 - is in princi- 
ple able to distinguish between a Kolmogorov-like or a different turbulent regime, 
determining 5 with an uncertainty of about 10-15%, if systematic errors are low 
enough. If the latter were important and at a level of about 10%, the uncertainties 
on the derived S would be huge. Experiments reaching collecting powers of hundreds 
m 2 sr days have the potential of fixing 5 at about 85% confidence level even in the 
case of important systematic errors, or at > 90% if the latter are of the order of the 
statistical ones. 

Only space-based or satellite detectors will be able to determine 5 with very high 
accuracy. In fact, even if they would be affected by not negligible systematic errors, 
the covered energy range would be significantly wider (we discussed the case with 
maximal energy of order 30 TeV/n). 

A careful measurement of the single fluxes will shed new light on the acceleration 
spectrum: in the high energy region the new generation apparata will be able to 
measure the spectrum slope 7 for light nuclei with a negligible error. No further un- 
certainty other than the one affecting 5 as derived from the B/C ratio will then be 
added to the determination of the acceleration slope a. This last parameter is thus 
expected to be measured with an uncertainty of 10-15% in the near future, allowing 
to exclude large ranges of a values and leading to a much clearer understanding of 
the acceleration processes. 

The spectra of the single elements of the CR beam will be measured with very good 
statistics up to more than some 10 14 eV (which is the correct unit to be used when 
comparing with indirect measurements, for which only the energy /particle is deter- 
mined). With a careful analysis of the experimental systematics, absolute particle 
fluxes will be provided, thus giving more than a decade calibration region to the 
Extensive Air Shower detectors. Direct measurements based on adequate statistics 
in the region of overlap with the Extensive Air Shower explored region are of main 
importance in order to give ground based experiments a firm reference point, thus 
helping in deriving the mean composition of the CR beam towards and at the knee 
and in checking the models which are used to derive the energy. Around 10 14 eV, the 
flux uncertainty can be set at present around 30%; with the new experiments now 

4 http:/ /cosmicray.umd.edu/cream/cream.html 
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under way or in project, we can expect to significantly reduce it. The same is true 
for what regards the mean mass: results from indirect measurements are far from 
clear [51], while data from direct detection still show about 30% uncertainty above 
10 14 eV, which can surely be lowered at least to 15-20% with the new apparata. 
At high energy, above 100 GeV/n, the primary nucleons contribute mainly to 
neutrino-induced upward through-going or stopping muons [52]; most important 
are protons and helium nuclei, while heavier ones contribute less than ~10% at all 
energies. Apparata with the considered exposures will be able to measure also pro- 
ton and helium spectra with unprecedented accuracy, thus allowing to reduce the 
present uncertainty in the all nucleon spectrum, which including all the available 
data is still around 30% in the TeV range. 

We want to emphasize that a better understanding of the propagation properties 
is of great significance for the study of signatures of new physics in CRs. Indi- 
rect signals of dark matter pairs annihilating in the halo of our Galaxy could be 
found in antiproton, antideuteron or positron CRs. This research is limited by the 
uncertainties in the propagation parameters, that dramatically affect the fluxes of 
charged particles located in the whole diffusive halo [53]. For instance, if this uncer- 
tainty were significantly reduced, it would be possible to exclude supersymmetric 
models predicting antiproton fluxes exceeding the present measurements [54]. A 
better signal-background discrimination would come from antideuteron measure- 
ments [55], and a reduction of the astrophysical uncertainties on the calculated 
fluxes would be even more desirable. 

We conclude by reminding that important parameters - other than 5 and a - de- 
scribing our Galaxy in relation to the propagation of charged particles can be fixed 
only with very precise low energies data, and for various species. This in princi- 
ple could be done by the AMS experiment on the Space Station [26], which could 
measure fluxes from hundreds MeV/n up to TeV/n and for 1< Z < 26. 
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